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Neutrino properties can be constrained by the detection of ultra-high energy cosmic neutrinos 
(UHECNs). By using the updated global fitting results of neutrino mixing parameters, we present 
predictions on the neutrino flavor ratios at the Earth from three possibly astrophysical sources. 

Comparing with the latest IceCube data, we find that the normal hierarchy (NH) and inverted 
hierarchy (IH) cases from the initial ratios </>° e : q^ =1:2:0 and 0:1:0 are compatible with 

the data in the standard neutrino oscillation scenario. We also examine the neutrino flavor ratios 
in a neutrino decay scenario beyond the standard model, and introduce the special case that two 
mass eigenstates of neutrinos, i.e., and v%, are degenerated. We find that the IH case and the 
degenerate NH case from the 1:2:0 and 0:1:0 sources are still permissible with the IceCube data 
within the 3cr error range. The general constraints only rely on the neutrino mixing and oscillation 
framework are also discussed. 
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The properties of neutrinos are under active investiga¬ 
tions from various experiments. The neutrino oscillation 
has been well established by the experiments with solar, 
atmospheric, reactor and accelerator neutrinos in recent 
decades [1-9]. There have been a number of new mea¬ 
surements on the neutrino mass splitting recently [10-13]. 
However, the hierarchy of mass eigenstates is still unde¬ 
cided, and it is also not clarified yet whether neutrinos 
are Dirac or Majorana fermions. 


The detection of ultra-high energy cosmic neutrinos 
(UHECNs) can provide us information concerning neu¬ 
trino properties [14, 15]. There are many running or on¬ 
going experiments such as ANITA [16], ARA [17], and 
IceCube [18]. The IceCube collaboration just reported 
the latest results of ultra-high energy cosmic neutrino 
flavor ratios [19, 20]. The UHECNs with energy above 
10 7 eV are expected to arise from a series of processes 
in the ultra-high energy cosmic rays (UHECRs). One of 
the components of the UHECRs is the ultra-high energy 
proton flux, which is thought to origin from the extra- 
galactic processes such as gamma-ray bursts (GRBs) or 
active galactic nuclei (AGNs) [21]. These accelerated pro¬ 
tons may scatter with the intergalactic medium or the 
cosmic microwave background (CMB) photons and gen¬ 
erate the pions: 

p+p —> n + TT~n °, (1) 

p + 7 -A A + —> n + tt + . (2) 

The proton-photon scattering in Eq. (2) is known as the 
Greisen-Zatasepin-Kuzmin (GZK) process [22, 23]. The 
charged pions and their muon daughters generate the 
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UHECNs by the decay chains 

7 r + —> —> e + + Vp + 9^ + v e , (3) 

7T - -t H~ + P/j, -t e~ + + v e , (4) 

which result in the most conventional initial ratios in the 
pion-decay case [24, 25] 

< : < : = 1 : 2 : 0. (5) 

Besides the pion decays, some other generative mecha¬ 
nisms of UHECNs were proposed. If the muons lost en¬ 
ergy while they are going through strong magnetic fields 
or matters [26, 27], the energy of is e and v e emitted in the 
secondary decays in Eqs. (3) and (4) are lower than the 
ultra-high energy range we concern. Hence, at high en¬ 
ergy the v e flux is decreased and the ratios approximately 
change into the muon-damped case: 

= 0 : 1 : 0 . ( 6 ) 

On the other hand, the decays of neutrons from the heavy 
nuclei photo-dissociation provide a pure v e generation to 
bring the ratios in the neutron-beam case [28] : 

< :<.:<= 1:0:0. (7) 

The flavor ratio detection can provide us information con¬ 
cerning the UHECN sources. 


Since the neutrinos are mixing among the mass eigen¬ 
states in the flavor eigenstates, the neutrino flavor ra¬ 
tios change during propagating and the final ratios at 
Earth are different from the initial ratios at sources. 
In the three-generation neutrino framework, the neu¬ 
trino mixing is well described by a 3 x 3 unitary matrix 
U which is known as the Pontecorvo-Maki-Nakawaga- 
Sakata (PMNS) mixing matrix [29, 30]. The flavor mix- 
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ing of three neutrino generations is written as 

Wi)=J2 u ijWj), ( 8 ) 


states and j = 1, 2, 3 represent the mass eigenstates. In 
the standard parametrization, the PMNS matrix is ex¬ 
pressed by three mixing angles 012,013,023 and one CP- 
phase angle <5 in a form 


where the subscripts l = e, /x, r denote the flavour eigen- 


V{6\2i $ 23 , ^ 13 , 8) 


C12C13 

— Ci2S23Sl3e I<5 — S12C23 

— Ci2C23Si3e 1<5 + S12S23 


S12C13 

— Si2S23Sl3e 1<5 + C12C23 

— Si2C23Si3e 1<5 — C12S23 


si 3 e- w \ 
S 23 C 13 I , 
C 23 C 13 / 


(9) 


TABLE I: The global fit of neutrino mixing parameters [31], 
with NH or IH denoting the normal or inverted hierarchy of 
mass eigenstates. 


parameter best fitiler 3er range 


sin 2 0 i 2 (NH or IH) 0.304±g;gi| 
sin 2 6*23 (NH) 0.452±g;g|| 

sin 2 0 23 (IH) 0.579^0^25 

sin 2 0 13 (NH) 0.0218to'ooio 

sin 2 0 13 (IH) 0.0219to'oon 

5/°(NH) 306+70 

<y/°(IH) 254+02 


0.270 -A 0.344 
0.382 -> 0.643 
0.389 0.644 

0.0186 ->• 0.0250 
0.0188 ->• 0.0251 
0 ->• 360 
0 ->• 360 


where s l3 and c i3 denote sin0 i:) and cos0^ ( i,j = 1,2,3). 
Here we use the recently updated global fit data [31] 
listed in Table I. 

For an individual neutrino arisen at a determined po¬ 
sition, the transition probability is related to the mass 
square difference and the propagation length. But if per¬ 
forming an average over all the neutrinos arisen at ar¬ 
bitrary positions, only mixing parameters are related to 
the ratio changes. The flavor ratio detected at Earth q i® 
is written as [25] 


^ = Ei^i 2 i^i 2 ^- ( 1Q ) 

V ,i 

If considering a tribimaximal mixing case, the flavor ra¬ 
tios in Eq. (5) originated from the pion-decay sources 
change into [25] 


= ( n ) 

Using the global fitting values of mixing parameters, we 
can evaluate the final detected flavor ratios of different 
initial ratios discussed above. The results are shown in 
Table II and Fig. 1. The ratios in a determinate mass 
hierarchy are in a line, because the ratios at source are 


TABLE II: The predicted flavor ratios of neutrinos at the 
Earth in different cases and initial ratios. The ratios are in 
the v e : 1 /^ : v T order. The best-fit data of the IceCube 
experiment can be found in Refs. [19, 20] 


cases ratios at source ratios at the Earth 



1 

2 

0 

0.35 : 

0.33 

0.32 

(NH) 





0.31 : 

0.35 

0.34 

(IH) 

standard 

0 

1 

0 

0.25 : 

0.36 

0.39 

(NH) 

oscillation 




0.19 : 

0.43 

0.38 

(IH) 


1 

0 

0 

0.55 : 

0.25 

0.20 

(NH) 





0.55 : 

0.19 

0.26 

(IH) 

decay (NH) 




0.68 : 

0.21 

0.11 


decay (IH) 




0.02 : 

0.57 

0.41 



1 

2 

0 

0.50 : 

0.27 

0.23 


deacy 

0 

1 

0 

0.47 : 

0.28 

0.25 


(degenerate) 

1 

0 

0 

0.56 : 

0.25 

0.19 


IceCube 




0.49 : 

0.51 

0 (new) 

best-fit 




0 : 0 . 

2 : 0.8 

(old) 



linear correlated, i.e., 

(1 : 2 : 0 ) = (1 : 0 : 0)+2 x (0 : 1 : 0 ). ( 12 ) 

The ranges marked by the solid lines in Fig. 1 represent 
the error ranges due to the 3cr range of the mixing param¬ 
eters. Comparing with the latest data from the IceCube 
experiment [19, 20], we find that in both NH and IH 
cases, the pion-decay sources and muon-damped sources 
are compatible with the IceCube data, but the neutron- 
beam sources are disfavored even thought we broaden 
the error range to 3cr. Since the flavor ratios in different 
mass hierarchy cases have distinct deviations, it is ex¬ 
pected that the future data can put strong constraint on 
the mass hierarchy in the standard oscillation scenario. 

Now we consider a new scenario beyond the standard 
model. Since the splitting among three mass eigenstates 
have been detected, the heavier neutrinos may decay into 
the lighter ones via flavor changing processes [32-34], 
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propagation, i.e., L r, only the lightest neutrino can 
arrive at Earth. Thus, the final flavor ratios detected 
at Earth depend on the mass hierarchy and the mixing 
parameters. The mass hierarchy determines the remain 
mass eigenstate, and the mixing parameters determine 
the flavor components of each mass eigenstate. If the 
neutrino mass eigenstates are in normal hierarchy (NH), 
i.e., m 3 > m 2 > mi, the UHECN beam contains the only 
stable mass state zq. So that Eq. (10) reduces to 


4>f = \Un\ 2 . (14) 


The flavor ratios at Earth are 



= 0 . 681 “;“ 


0.21 


+ 0.08 

- 0.16 


0.11 


+ 0.17 

- 0 . 06 - 


(15) 


FIG. 1: The comparison of the predicted flavor ratios at Earth 
in the standard oscillation scenario with the IceCube data. 
The NH cases are marked by the solid points corresponding to 
the initial ratios by shapes, and the IH cases are marked by the 
open points. The regions marked by the solid lines represent 
the error ranges due to the 3 ct range of mixing parameters. 
The 68% and 95% confidence regions of the IceCube data are 
also indicated. 
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FIG. 2: The comparison of the predicted flavor ratios in the 
the neutrino decay scenario with the IceCube data. The non¬ 
degenerated cases in NH and IH are marked by the inverted 
triangles. The points of the degenerated cases from different 
sources are in a line, for the similar reason of the standard 
oscillation scenario in Fig. 1. 


The dominate neutrino decay manner is thought to be 

zq — y Vj T A and v j — y Vj T A., (13) 

where zq [i = 1,2,3) are neutrino mass eigenstates 
and A is a very light or massless particle, e.g., a Ma- 
joron [14], Many more detailed discussions have been 
proposed [15, 35-38]. Since the exiting limit on the half- 
life r of the decay channel in Eq. (13) is too weak to elim¬ 
inate the UHECN decay with high energy and lengthy 


The 3cr range of the mixing parameters is considered in 
the error calculation. A probable ratio is an arbitrary but 
unitary combination of the three fa's in the error range, 
as shown in Fig. 1. If the neutrino mass eigenstates are 
inverted hierarchy (IH), i.e., m 2 > mi > m 3 , the remain 
mass eigenstate turns to be zq. Hence the flavor ratios 
at Earth change into 



€ t = 0.02 


+ 0.003 

- 0.003 


0.57 


+ 0.06 

- 0.19 


0.41 


+ 0.019 
- 0.06 • 


(16) 

Further more, we consider the neutrino decay in de¬ 
generate cases. Since the mass splitting between zq and 
zq, i.e., Arn^i, is two orders of magnitude smaller than 
Am§! and Am § 2 [ 10 - 12 ], it is possible to conjecture that 
the mass eigenstates zq and zq are approximatively de¬ 
generate relative to zq [39]. This assumption is supported 
by the experiments of the neutrino-less double beta de¬ 
cay [39] and the cosmological observations [40] , since the 
mass splitting Am]j is much smaller than the bounds of 
absolute mass scales. For the inverted hierarchy, both zq 
and zq decay to zq, so still only the zq is stable. Hence, 
the degenerate IH case leads to the same result as that 
of the non-degenerate IH case. But in the normal hier¬ 
archy case, zq is the heaviest mass eigenstate, so both zq 
and zq are stable due to the degeneration. Each stable 
mass eigenstate originates from not only the initial flavor 
eigenstates but also the zq decay. Since the flavor com¬ 
ponents of zq and v 2 are disparate, i.e., Un 7 ^ Ui 2 1 and 
the branching ratios of the zq decay may be different, the 
latter two factors in Eq. (10) can not be cancelled. That 
is to say, the final flavor ratios at the Earth are related 
to the initial flavor ratio source and the branching ratios. 
As a straightforward assumption, i.e., letting the branch¬ 
ing ratios of zq and zq both being equal to 50%, Eq. (10) 
changes to 


= E \Uu\ 2 (\Uih\ 2 + \Uu 3 \ 2 /2)4, (17) 

V ,i 

where i = 1,2 and V = e, //, r. Then we can get different 
final flavor ratios for different sources. 

The results of neutrino decay scenarios are shown in 
Table II and Fig. 2. The 3a error regions of degenerate 
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NH cases are smaller than those of non-degenerate cases, 
because the flavor ratios in Eq. (17) are more complicated 
and the errors from different mixing parameters cancel 
each other, due to the unitary of the PMNS matrix. The 
error region of the IH case is narrower along the v e di¬ 
rection, because the error range of $13 is much smaller 
than those of other mixing parameters. We can find that 
the IH case (only 1/3 is stable) and the degenerate NH 
case from the pion-decay sources ( 1 : 2 : 0 ) and the muon- 
damped sources (0:1:0) are compatible with the IceCube 
data at the 3cr error range, although the best-fit points 
in these cases are at the boundary of 95% confidence re¬ 
gion. The non-degenerate NH case and the degenerate 
NH case from the neutron-beam sources (1:0:0) are dis¬ 
favored comparing with the latest IceCube data. 

In the early years, the flavor ratios of UHECNs were 
discussed on a premise that the matrix element t /13 
is vanishing [14, 15]. After the discovery of the non¬ 
vanishing mixing angle $13 [9], many theoretical calcula¬ 
tions based on the updated mixing parameters are pro¬ 
posed [41-49], in both standard and non-standard sce¬ 
narios. According to the new result announced by the 
IceCube collaboration [20], two kinds of sources in the 
standard oscillation scenario, i.e., pion-decay sources and 
muon-damped sources, are well compatible with the new 
IceCube data, while the non-standard oscillation scenar¬ 
ios have not been found. We propose a new probable 
case that the mass eigenstates mi and m 2 are degener¬ 
ate in the neutrino decay scenario. The degeneracy of 
mass eigenstates changes the stable states, and results in 
different final flavor ratios from different sources. In the 
degenerate NH case, we find that the theoretical predic¬ 
tions of the neutrino decay scenario are still permissible 
comparing with the IceCube data, by using the updated 
global fitting mixing parameters. 
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FIG. 3: The comparison of the predicted flavor ratios in the 
arbitrary initial ratio scenario with the IceCube data. The 
best fit regions in NH and IH are the line shape region and 
the thin triangle in the middle area. The 3er region is denoted 
by the filled area. The IceCube best fit points, 68% and 95% 
confidence regions are also shown as comparisons. 


In the discussions above, the initial flavor ratio at the 
source was introduced as a constraint. Different initial 
ratios lead to different final ratios that are detectable 
by the IceCube experiment. Hence, one can conjecture 
which source or sources play the major roles according to 
the IceCube data or which way of the neutrino decay is 
probable, according to the IceCube data. On the other 
hand, if all the possible initial ratios are included, i.e., 
the initial ratios are arbitrary, one can also get a con¬ 
straint on the final ratios by the restrict on the neutrino 
mixing and oscillation only. Here we consider a general 
example that not only flavor ratio is arbitrary, even the 
ratio of mass eigenstates is also arbitrary, which means 
the mass eigenstates are allowed to be incoherent during 
propagating. The later case is more general since it can 
completely cover the region of the former case. Fig. 3 
shows the final ratios with the mixing parameters best 
fit regions and the 3a regions in both NH and IH case, 
compared with the IceCube data. The 3cr regions in NH 
and IH case are almost the same. In Fig. 3, the best fit 
regions makes up only a very small portion in the ratio 
triangle, no matter in NH or IH case. It suggests that 
the neutrino mixing framework can actually provide a 
strong constraint on the flavor ratio detected at Earth, 
no matter what happened at the source. 

Another point can not be neglected is that the best fit 
points of the IceCube data are not covered by the general 
3cr region, though the 95% and 68 % confidence regions 
overlap with the 3er region in some areas. That may 
indicate the probability of new physics models beyond 
the the neutrino mixing and oscillation framework today. 
It is well expected that the UHECNs detections can be 
used as constraints of the mixing parameters and the 
mechanism. 

In summary, we discuss the flavor ratios of UHECNs 
originated from three possibly astrophysical sources. The 
final flavor ratios at the Earth depend on the initial ratios 
because of the neutrino oscillation or the neutrino de¬ 
cay. By using the recently updated global fitting results 
of mixing parameters, we evaluate the final flavor ratios 
within the 3er error range, and compare our results with 
the recent IceCube data [19, 20]. In the standard oscilla¬ 
tion scenario, both NH and IH cases from the pion-decay 
sources ( 1 : 2 : 0 ) and the muon-damped sources ( 0 : 1 : 0 ) are 
compatible with the IceCube data, but the neutron-beam 
sources (1:0:0) are disfavored for both NH and IH cases. 
If considering the neutrino decay, we find that the IH 
case and the degenerate NH case from the 1:2:0 and 0:1:0 
sources are permissible within the 3er error range, while 
the non-degenerate NH case or the degenerate NH case 
from the 1:0:0 sources are disfavored with the data. If 
removing all the other restricts but neutrino mixing and 
oscillation, we can also get a constraint on the final ra¬ 
tios. Since the difference between NH and IH cases are 
considerably significant in the standard scenario, the de¬ 
cay scenario and the general constraint, the future data 
of UHECNs can provide us more information on the neu¬ 
trino properties. 
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